BARI-TH 453/03 



Dynamics of Ferromagnetic Walls: Gravitational Properties 



L. Campanelli 1 ' 2 *, P. Cea 3 ' 4t G. L. Fogli 3 ' 4 * and L. Tedesco 3 ' 4§ 

1 Dipartimento di Fisica, Universita di Ferrara, 1-44100 Ferrara, Italy 
2 INFN - Sezione di Ferrara, 1-44100 Ferrara, Italy 
3 Dipartimento di Fisica, Universita di Bari, 1-70126 Bari, Italy 
4 INFN - Sezione di Ban, 1-70126 Ban, Italy 



Abstract 

We discuss a new mechanism which allows domain walls produced during the 
primordial electroweak phase transition. We show that the effective surface tension 
of these domain walls can be made vanishingly small due to a peculiar magnetic 
condensation induced by fermion zero modes localized on the wall. We find that 
in the perfect gas approximation the domain wall network behaves like a radiation 
gas. We consider the recent high-red shift supernova data and we find that the 
corresponding Hubble diagram is compatible with the presence in the Universe of a 
ideal gas of ferromagnetic domain walls. We show that our domain wall gas induces 
a completely negligible contribution to the large-scale anisotropy of the microwave 
background radiation. 
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1. Introduction 



A considerable amount of interest has emerged in the physics of topological defects pro- 
duced during cosmological phase transitions [T] . Indeed, even in a perfectly homogeneous 
continuous phase transition, defects will form if the transition proceeds sufficiently faster 
than the relaxation time of the order parameter [21 El IU IHj ■ In such a non-equilibrium 
transition, the low temperature phase starts to form, due to quantum fluctuations, simul- 
taneously and independently in many parts of the system. Subsequently, these regions 
grow together to form the broken-symmetry phase. When different causally disconnected 
regions meet, the order parameter does not generally match and a domain structure is 
formed. 

In the standard electroweak phase transition the neutral Higgs field is the order parameter 
which is expected to undergo a continuum phase transition or a crossover [HIE!- In the 
case in which the phase transition is induced by the Higgs sector of the Standard Model, 
the defects are domain walls across which the field flips from one minimum to the other. 
The defect density is then related to the domain size and the dynamics of the domain 
walls is governed by the surface tension o. The existence of the domain walls, however, 
is still questionable. It was pointed out by Zel'dovich, Kobazarev and Okun [H] that the 
gravitational effects of just one such wall stretched across the universe would introduce 
a large anisotropy into the relic blackbody radiation. For this reason the existence of 
such walls was excluded. Quite recently, however, it has been suggested EH that the 
effective surface tension of the domain walls can be made vanishingly small due to a pe- 
culiar magnetic condensation induced by fermion zero modes localized on the wall. As 
a consequence, the domain wall acquires a non zero magnetic field perpendicular to the 
wall, and it becomes almost invisible as far as gravitational effects are concerned. In a 
similar way, even for the bubble walls, which are relevant in the case of first order phase 
transitions, it has been suggested [TT] that strong magnetic fields may be produced as a 
consequence of non vanishing spatial gradients of the classical value of the Higgs field. 
It is worthwhile to stress that in the realistic case where the domain wall interacts with 
the plasma, the magnetic field penetrates into the plasma over a distance of the order of 
the penetration length, which at the epoch of the electroweak phase transition is about 
an order of magnitude greater than the wall thickness. This means that fermions which 
scatter on the wall feel an almost constant magnetic field over a spatial region much 
greater than the wall thickness. 

In this paper we consider thin domain walls eventually produced during the primordial 
electroweak phase transition. We focus on domain walls where the magnetic condensa- 
tion induced by fermion zero modes leads to ferromagnetic domain walls (in the following 
FDW). Obviously, ferromagnetic domain walls are not topological stable. However, FDWs 
are characterized by the vanishing of the effective surface tension. As a consequence FDWs 
are dynamically stable structures, due to a huge energy barrier against spontaneous de- 
cays, which could have survived until today. 

The plan of the paper is as follows. In Section II we consider FDW in the thin wall approx- 
imation and discuss the general properties of the associated energy-momentum tensor. In 
Section III we evaluate the energy momentum tensor contributions due to positive energy 
fermions captured by the domain wall. In Section IV we discuss the equation of state 
of a domain-wall network in the ideal gas approximation. In Section V we discuss the 
gravitational properties of the ideal gas of FDWs. Section VI is devoted to the analysis 
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of the temperature anisotropy on the Cosmic Microwave Background Radiation induced 
by the ferromagnetic domain wall network. Finally, we summarize our results in Section 
VII. Some technical details are relegated in the Appendix. 

2. The energy-momentum tensor of a FDW 

Topologically stable kinks are ensured when the vacuum manifold A4, the space of all 
accessible vacua of the theory, is disconnected, i.e. n (M.) is nontrivial. 
We shall consider a simplified model in which the kink is a infinitely static domain wall in 
the (x, y) plane in a flat space-time. That is we assume that the vacuum manifold consists 
of just two disconnected components and restrict ourselves to one spatial dimension. In 
our model the scalar sector that give rise to a planar wall is given by a real scalar field 
with density Lagrangian, 

c^= l -d^d^ -V{\<j>\,T). (2.1) 

Lagrangian (2.1) has an exact Z2 symmetry corresponding to the discrete transformation 
(f) — > —(f). At zero temperature the Z2 symmetry is spontaneously broken, with <fi acquiring 
a vacuum expectation of order v. In the tree approximation, the set of vacuum states is 
((f)) 2 = v 2 so that n (M) 7^ 0. Then, if we apply boundary conditions that the vacuum 
at —00 lie in — v and that at +00 lie in +v, by continuity there must exists a region in 
which the scalar field is out of the vacuum. This region is a domain wall. If we take the 
potential at zero temperature as V((f>) = j((fi 2 — v 2 ) 2 , it is easy to see that the classical 
equation of motion admits the solution describing the transition layer between two regions 
with different values of (0), 

(f){z) = v tanh(z/A), (2.2) 

where A = \f2\~ 1 ^ 2 v~ l is the thickness of the wall. 
The energy-momentum tensor for the 0-kink (2.2) is, 

T^ )v = a{z)a^diag{l,l,\^), (2-3) 

a(*0 = ^[cosh(z/A)]- 4 , (2.4) 

where 

a, = |^T ( V = ^/V (2 . 5) 

is the surface energy density associated with the scalar field. In the thin wall approxi- 
mation, i.e. A — > 0, a(z) reduces to a 5-function. Unless the symmetry breaking scale 
v is very small, surface density energy of the kink is extremely large and implies that 
cosmological domain walls would have an enormous impact on the homogeneity of the 
Universe. A stringent constraint on the wall tension for a Z2-wall can be derived from 
the isotropy of the microwave background. If the interaction of walls with matter is negli- 
gible, then there will be a few walls stretching across the present horizon. They introduce 
a fluctuation in the temperature of the microwave background of order 5T/T ~ G^a^to, 
where to is the present time. Observations constrain 5T/T < 1CT 5 , and thus models 
predicting topologically stable domain walls with v > 1 MeV should be ruled out. 
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On the other hand, we are interested in domain walls formed during the electroweak phase 
transition by the Kibble-Zurek mechanism [21 El EE]- Even these domain walls are ap- 
proximately accounted for by Eq. (2.2), so that the energy-momentum tensor is given by 
Eq. (2.3). Since v ~ 10 2 GeV, it is widely believed that there must be some provision for 
the elimination of these defects which indeed are not topological stable. However, quite 
recently it has been suggested that the effective surface tension of the domain wall can 
be made vanishingly small due to a peculiar magnetic condensation induced by fermion 
zero modes localized on the wall. As a consequence the domain wall acquires a non zero 
magnetic field perpendicular to the wall. This mechanism is able to make these defects 
almost dynamical stable. Indeed, FDWs are protected against decay by a huge energy 
barrier. 

In general, for a FWD the magnetic field vanishes in the regions where the scalar conden- 
sate is constant. So that, the magnetic field can be different from zero only in the regions 
where the scalar condensate varies, i.e. in a region of the order of the wall thickness. How- 
ever, it is worthwhile to stress that in the realistic case where the domain wall interacts 
with the plasma, the magnetic field penetrates into the plasma over a distance of the order 
of the penetration length which at the time of electroweak phase transition is about an 
order of magnitude greater than the wall thickness. This means that fermions interacting 
with the wall feel an almost constant magnetic field over a spatial region much greater 
than the wall thickness. So that as concern the interaction of fermions with FWDs, we 
can assume that the magnetic field is almost constant. 

For completeness, let us briefly review the generation of ferromagnetic domain walls at 
the electroweak phase transition j9]. After the electroweak spontaneous symmetry break- 
ing the unique long range gauge field is the electromagnetic field A^. In our scheme the 
real scalar field is the physical Higgs field, so that A^ does not couple directly to 0. The 
electromagnetic Lagrangian is: 

C em = - l -F flv F^. (2.6) 

We are interested in the case of an external magnetic field B (localized on the wall at 
z = 0) directed along the ^-direction, 

A M = (0,0,-Sx,0), (2.7) 

where without loss of generality we may consider B > 0. Assuming that the magnetic 
field B is localized on the wall, we can write, 

T ( B ) » = °B diag (1, -1, -1, 0), (2.8) 

where 

/+oo -q2 

dzT° B)0 = A— (2.9) 

is the surface energy density associated with the magnetic field if B is sizable over a 
distance z ~ A, and (3(z) is a localization function of B on the wall. In the thin wall 
approximation we expect that (3{z) reduces to a 5-function. 

Next we consider a four dimensional massless Dirac fermion coupled to the scalar field 
through the Yukawa coupling, and in presence of the homogeneous background magnetic 
field (2.7). The Dirac fermion if> is described by the Lagrangian: 

C^ = iP(ip-e4-g Y( t>)iP (2.10) 
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where we consider e > 0. In presence of the kink (2.2) the Dirac equation, 

(ip-e4-g Y cj ) )^ = 0, 



(2.11) 



admits solutions localized on the wall [Oj. Using the following representation of the Dirac 
algebra, 

' l a \ s ( % 



a = 0,1,2 : 7" 



- 7° 



7 



i 



with 



7° = <x 3 , 7 1 = icr u I 2 = ia 2} 



where a, are the usual Pauli matrices, we find for the localized modes [Oj, 



ip(x,y,z,t) = ^uj(z) ( P^y^) 



(2.12) 
(2.13) 

(2.14) 



lu(z) = N [ cosh (z/A)]^ A , (2.15) 

where £ = gy v is the mass which the fermion acquires in the broken phase, p(x, y, t) is a 
Pauli spinor that satisfies the (2 + l)-dimensional massless Dirac equation in a constant 
external magnetic field, and N is a normalization constant given by 



N 



dzto 2 (z) 



-1/2 



[AB(£A,l/2)] 



-1/2 



(2.16) 



B(x , y) being the Euler Beta function. 

Let us consider the energy-momentum tensor of the four dimensional massless Dirac 
fermion, T^ v . It is straightforward to show that the non-zero components are given by: 



T ( ^ = ^ 2 (z)pt7°7 Q (V, 



(2-17) 



% 3 = -igr <f>(^ 2 (z)p^ 7° 7> (2-18) 

In the thin-wall approximation uj 2 (z) reduces to a 5-function so that T^ 3 goes to zero. 
Therefore the tensor reduces to: 

Tfc v (x,y,z)=6(z) , (2.19) 



where Tf\p = ip^ 7° l a dpp is the energy-momentum tensor for a (2 + l)-dimensional 

massless Dirac fermion coupled to the magnetic field. 

The resulting energy-momentum tensor of the wall is given by 

Tfa v = %) v + Tfav + Tfav- (2-20) 

^From Eq. (2.19) we see that the (3 + l)-dimensional tensor v is obtained by evaluating 
which corresponds to (2 + l)-dimensional Dirac fermions localized on the wall. 



It is known since long time that in (2 + 1) -dimensional QED in presence of un external 
homogeneous background magnetic field with massless Dirac fermions it is energetically 
favorable the spontaneous generation of a negative mass term ^2]- So we are lead to 
consider the following effective (2 + l)-dimensional Dirac equation: 

(i Td a -eTA a -m) p(x, y, t) = 0. (2.21) 
Taking p(x, y, t) = e~ lEt p(x, y) we get [12] 

E = +K: W) \ (2.22) 



E = -E n : pS=('^V /2 ^e-^| . |, (2.23) 



2E n I J^k 



N n H n (C) 



E n +m 



E = E« = m : /;,,„ = A e~^ 2 ( J ). (2.21.) 



where 



E n = V2neB + m 2 , (2.25) 
C = sfeB (x - -|) , (2.26) 

N , = (^^, ( , 27) 



V 7T / y/2'-n\ 

with n = 1,2, and H n (x) are Hermite polynomials. The normalization is 



/ ^pg^x)^ (x) = 6{p-p')6 nn ,. (2.28) 



It is convenient to expand the fermion field operator p in terms of pn& and pi, J . In other 
words, we adopt the so-called Furry's representation. For negative mass term we have: 



P 



E / + E / d P b npPn} ■ ( 2 - 29 ) 



We observe that in the case of negative mass the positive solutions have eigenvalues +E n 
with n > 1 and the negative ones — E n with n > 0. In the expansion of p we associate par- 
ticle operators to positive energy solutions and antiparticle operators to negative energy 
solutions. These operators satisfy the standard anticommutation relations, 

i a n P ,4'p'l = i b np ,bl p >} = $nn>$(p-p'), (2.30) 

others anticommutators being zero. The energy-momentum tensor operator is 

r ( J )/J = ip t 7°7 a ^p. (2.31) 
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We need the surface density of the vacuum expectation value of energy-momentum tensor 
operator defined as: 

^10) 



(0|0) 



(2.32) 



where A is the area of the wall and |0) is the fermion vacuum. In the Furry's representation 
(T$ )fl ) reads 

< T (p)/3> = T /^E / dppWiVl-dpptt. (2.33) 

J n=0 J -°° 

for the negative mass theory. Inserting Eqs. (2.23) and (2.24) into Eq. (2.33), it is 
straightforward to obtain 



(p)o 



eB °° 



n=0 



P P = (^)i) = (T 



p 2 , 



(eB) sr-^ n 
2tt ^ K ' 

n=0 " 



(2.34) 
(2.35) 



while all other ( ^ ) are equal to zero. Therefore, the energy-momentum tensor is: 



( T (p) a ) = dia 9 (°> ~Pp, -P P ), 



(2.36) 



where a p and p p are the surface density energy and pressure of the fermionic condensate, 
respectively. It easy to show that the following relation 



P P = B 



d(Jp 



holds. We shall consider the "effective" energy-momentum tensor defined as: 

( T (p) 3) 'ef f = (T(p)f3)B -{T( p )p)b=0- 

By using the integral representation 

r0 ° ds d „. 



TTS ds 



introducing the dimensionless variable r\ = eB/m 2 and the function 



(2.37) 



(2.38) 



^2.39) 



ds d 
its ds 



-sir] 



1-e 



-2s 



2s 



(2.40) 



we cast Eq. (2.38) into 



= (TLo)eff = S(Z) 



(eB) 



3/2 



-Pf - (Tin) 



(p) 1 /eff 



2tt 

(p)2)eff 



g(v), 



(TL 2 )eff = -s(z) 



(eB) 3 / 2 



2tt 



1 / x . d 9(v) 



(2.41) 
(2.42) 
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Taking into account Eqs. (2.20), (2.41), (2.42), the final expression of the effective energy- 
momentum tensor of a thin FDW is: 

( T (w) v ) = S ( z ) dia 9 ( 0v> , ~Pw , ~Pw , 0) , (2.43) 

where 

a w = a^ + a B + cr e / f , (2.44) 

p w = -a^ + a B +p e / f . (2.45) 

We are interested on the value of (T^ u ) for the magnetic field B* that minimizes the 
energy density a w . As shown in Ref. [5| the magnetic field B* can be approximated by: 



B* = 





m 


2 




1 


47T 


771 





(2.46) 



Moreover the stability condition 

*■* = (2.47) 

gives in turns the value of \m*\. Equation (2.47) assures that our ferromagnetic domain 
walls have a vanishing effective surface tension. Moreover, it is remarkable that Eq. (2.47) 
implies: 

p w (B*,\m\*) = 0. (2.48) 

Note that Eq. (2.48) is an exact consequence of the stability condition Eq. (2.47) and 
does not rely on the approximate expression for B* given by Eq. (2.46). By using the 
standard values v ~ 250 GeV, oiqed — e 2 /47r ~ 1/137 and assuming A ~ 0.5 we obtain 
from Eqs. (2.46) and (2.47) the following value of the magnetic field [§], 

B* ~ 5.0 10 24 Gauss. (2.49) 

Interestingly enough, such a value of the magnetic field at the primordial electroweak 
phase transition could be relevant for the generation of the primordial magnetic field JTBI . 
It is worthwhile to stress that Eqs. (2.47) and (2.48) imply that our ferromagnetic do- 
main walls are almost invisible as concern the gravitational effects. However, as already 
discussed, in the realistic case where the domain wall interacts with the plasma, the 
magnetic field penetrates into the plasma over a distance of the order of the penetration 
length A which is about an order of magnitude greater than A. It follows that the above 
estimate of the induced magnetic field at the electroweak scale B* is reduced by a factor 

~ 0.3 which is still of cosmological interest [14]. Moreover, on the ferromagnetic 
domain walls there are also positive energy states. As a consequence fermions incident 
on the wall with an energy equal to the empty states can be trapped on the wall giving 
a non trivial contribution to the energy-momentum tensor. Thus, in order to investigate 
the gravitational properties of the ferromagnetic domain walls it is important to take care 
of these contributions to the energy-momentum tensor. 

3. The energy-momentum tensor of a massive FDW 

On the wall there are also available positive energy bounded states. Therefore when 
incident fermions have the same energy as the allowed bounded states on the wall, they 
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could be captured. In this case, the fermions bounded to the ferromagnetic domain wall 
do contribute to the energy-momentum tensor of the wall. In the following we shall refer 
to FDW bounded with positive-energy fermions as massive FDW. In this Section we 
calculate the energy-momentum tensor of massive FWDs. 
To discuss the localized solutions with positive energy we rewrite Eq. (2.11) as: 

(ip-e4-tg)ij)(x,y,z,t) = 0, (3.1) 

where £ = gyv, g(z) = tanh(^) and ip(x,y, z,t) is the wave function of the fermion 
eventually captured by the FDW and having positive energy corresponding to a allowed 
bound state on the wall. For reader convenience we relegate in the Appendix the details of 
our calculations. As shown in the Appendix, the wave functions of the localized fermionic 
modes are given by: 

= NA n e-^ +ip yy- iE \cosh^ A [H n Eu 1 + + 2n^B\H n ^u 2 + }, (3.2) 

tP 2 +e = N A n _ l e-^ +l vyy- lE \ C osh^y^{ Hn _ 1 Eul + ^^Hnull (3.3) 

i>\ = NA^e-^y-^^cosh^y^lH^Eul-^BlH^l], (3.4) 

= NA n e-^ +i ^- iE \cosh^ A [H n Eul-2n^/\^\H n ^ 1 u 1 + ], (3.5) 

where N is the normalization constant evaluated in the Appendix, and 

E n = ^2n\e\B, n=l,2,... (3.6) 

The wave functions (3.2), (3.3) correspond to fermions, while (3.4), (3.5) to antifermions. 
Note that fermions have spin parallel to the magnetic field, while antifermions have spin 
antiparallel to the magnetic field. 

In a previous paper [T^] we investigated the scattering of fermions off walls in presence 
of a constant magnetic field In particular, we discussed fermion states corresponding 
to solutions of Dirac equation (2.11) which are localized on the domain wall. We found 
that there is total reflection for fermions with parallel and antiparallel spin at energies 
E\ — £ 2 = 2 n |e| B and E\ — £ 2 = 2 (n+ 1) |e| B, respectively. Note that the difference is 
due to the fermion mass term £ which, indeed, vanishes on the wall where the system is in 
the symmetric phase. ^From the physical point of view, we see that fermions with asymp- 
totically positive kinetic energy equal to a/2 n fej B for parallel spin, and a/2 (n+1) \e\ B 
for antiparallel spin, can be indeed trapped on the domain wall. As a consequence, we see 
that ferromagnetic domain walls immersed into the primordial plasma are able to capture 
fermions by filling the Landau levels up to n = N max . Obviously, N max depends on the 
plasma temperature T, for we must have E n < T. In our case, the temperature at the 
electroweak phase transition is of the order of 10 2 GeV. Taking into account Eqs. (2.49) 
and (3.6) we get N max ±= 1. 

Let us indicate with \w) the generic massive FDW, namely a ferromagnetic domain wall 
filled with fermions (+e) and antifermions (— e), with "spin up" (s = 1) and "spin down" 
(s = 2) respectively, up to Landau level N max . Thus, the expectation value of the 

^ Similar calculations have been presented in Ref. [T^] . 
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energy-momentum tensor on the state \w) is: 

(T" )- 1 /i 2 x (w|T( " )Jw) (3 7) 

VwJ-aJ** {w \ w) > ^ 

where 

r^ }l/ = i* 7 ^. (3.8) 

To evaluate (TfaJ) we expand the fermionic operator $ in terms of the states (3.2)-(3.5) 
as follows 

* = EE / ^[aK.^^V'+e'^+^K^,*)^-"'^]- (3-9) 

n a s 

In the thin wall approximation we find, 

(Tfc v (z)) = 8(z)a 9 diag(l, -1/2, -1/2,0), (3.10) 

where 

9 \p R*l 3 / 2 Armaa: 9 Ip R*l 3 / 2 

d^<2* W > = ^ii£-L_ ^ j£ v L - (3- 11 ) 

7T z — ' 7T 

71=1 



Taking into account the value of B* given in Eq. (2.49), we see that Eq. (3.11) implies 
that is comparable to given by Eq. (2.5), so that massive FDWs could display 
important gravitational effects which, eventually, could be in contrast with observations. 
Indeed, the cosmological electroweak phase transition should lead to a network of domain 
walls. However, solving for the cosmic evolution of a domain-wall network is quite in- 
volved. Nevertheless some essential features can studied by considering the dynamic of 
the domain- wall network in the ideal gas approximation [T7j . 



4. Ideal gas of FDWs 

In this Section we shall follow quite closely the analysis of Ref. to calculate the 
equation of state for a ideal gas of massive ferromagnetic domain walls. 
Let us consider a perfect gas of walls moving with velocity v in a box of volume V ^> £, 
where £ is a typical distance between walls. The perfect gas approximation amounts to 
neglect any dissipative effects due to the interaction of the walls. We first consider a planar 
wall in the x — y plane. Because of the symmetry, the energy-momentum tensor will be a 
function only of z. Following Ref. J7j we find that the average energy-momentum tensor 
of the wall gas is: 

where is given by Eqs. (3.10) and (3.11) for massive FDWs, and (L) is the average 
wall separation. If the walls are moving with average velocity v in the + z direction, we 
can obtain the energy-momentum tensor by boosting W^ u in ^-direction. Repeating the 
procedure for walls moving in the — z, ± x, ± y directions and averaging over the walls, 
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we see that the average energy-momentum tensor for an ideal gas of massive FDWs is 
diagonal. Moreover, we have: 

P» = (T™)) = gy7 2 , (4.2) 

p w = (T&) = (Tg) = = \p w , (4.3) 

where 7 = (1 — v 2 ) is the Lorentz factor and v is the mean velocity of a wall. In 
Eq. (4.2) a w , the surface energy density for a massive FDW, is given by Eq. (3.11). From 
Eq. (4.3) we see that the equation of state for a ideal gas of massive FDWs does not 
depend on the wall velocity and coincides with the equation of state of a radiation gas. 

It is clear that in the more realistic case in which the walls interact with the primordial 
plasma the equation of state for a gas of massive FDWs is expected to be more stiff than 
that of a radiation gas. Writing the pressure as p w = ap w we expect that a < 1/3 in the 
primordial Universe (say before the time of structure formation). On the other hand, at 
the present time the dissipative effects due to the interaction of the walls can be safely 
neglected and then the perfect gas approximation holds. 

In the following we shall work in the ideal case of non-interacting walls since we will 
concentrate on physical process occurred at relatively recent times (i.e. after the beginning 
of structure formation). 



5. Large scale gravitational properties 

In a recent paper we investigated the gravitational properties of thin planar massive fer- 
romagnetic domain walls in the weak field approximation jTHj. In order to understand 
the cosmological implications of FDWs we must consider the standard model of the Uni- 
verse which is based upon the Friedmann- Robert son- Walker cosmological model |17j . The 
Friedmann-Lemaitre equations for the expansions parameter R are: 

R 2 k 8nG N A 

#T + ^-^- Ptot + 3, (5-1) 



where R= dR/dt is the derivative respect to the cosmic time t, k is the curvature constant 
( k equals to —1, 0, +1 for a Universe which is respectively open, flat, and closed), and A 
is the cosmological constant. Equation (5.1) tells us that there are three competing terms 
which drive the expansion: a term of energy (matter and radiation), the cosmological 
constant, and a curvature term. Recent observations on temperature anisotropy of the 
cosmic microwave background ^JJ 1201 121] strongly indicate that the Universe geometry 
is very close to flat. So that, in the following we assume k = 0. We are interested in the 
matter-dominated Universe where the radiation energy density and the matter pressure 
are negligible. Allowing for a FDW gas with density p w and pressure p w related by the 
equation of state Eq. (4.3), the Friedmann-Lemaitre equations become: 

2 8nG N A 
H = -^—{pm + Pw) + -, (5.3) 
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2|+ H 2 = -SnG NPw + A, (5.4) 

where 

H = | (5.5) 
is the Hubble parameter. By introducing the following f2-functions: 

Vt m = — , (5.6) 

Pc 

n w = P - , (5.7) 

Pc 

= 3^2, (5-8) 

where p c = 3H 2 /8irG N is the critical density, we cast Eq. (5.3) in the form: 

fi m + fi w + fi A = 1. (5.9) 

By defining the following dimensionless variables: 

T=t/t , (5.10) 
r = i?/^o, (5.11) 

where the index refers to the present time, and using Eq. (4.3), the equation of state 
for the ideal gas of FDW, we write Eq. (5.4) in the form: 

= t H [Qmor- 1 + Q w0 r- 2 + ft A0 r 2 ] 1/2 . (5.12) 

In Figure 1 we compare the expansion parameter as a function of time for three different 
models. We have fixed Q m o = 0.3 and consider (Qao = 0.0 , Q w = 0.7) , (Qao = 0.7 , 
Q w0 = 0.0) and (fi^o = 0.6 , Q w q = 0.1). In the first case R behaves quite similar to the 
Einstein-de Sitter model. The standard cosmological model corresponds to Q m0 = 0.3, 
^ao = 0.7. From Fig. 1 we see that the three models are almost identical up to the 
present epoch ( t H < 1.7). To discriminate among models we compare the Friedman- 
Robertson- Walker magnitude-redshift relation with the combined low and high redshift 
supernova dataset. To this end, we recall that the luminosity-distance is: 

where z is the redshift relative to the present epoch, C is the intrinsic luminosity of the 
source, and T the observed flux. Setting dt{z) = HqDl{z), for Friedmann-Lemaitre 
models it is well known that: 



dL{z) = {1 + z) 7^ sin (^[wwJ ,ff! ' 0<0 ' (5 - 14) 

r z dz' 

dL (z) = (l + z)l WWQ if n M = 0, (5.15) 

d L (z) = (l + ^-^sinhfv 7 ^ / TTTKnA if ^o>0, (5.16) 
V"*o V Jo H \ z )/ H oJ 
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where the Hubble parameter is given by 

H(z) = H [Q m0 (l + z) 3 + Q A0 + Sl w0 {l + z) A ] 1/2 • (5.17) 
The apparent magnitude m of an object is related to its redshift by the following relation: 

m(z) = 5log w d L (z) +M, (5.18) 

where M. = M + 5 log 10 (Mpc/D^) + 25 is the Hubble-constant-free absolute magnitude, 
M is the absolute magnitude, and Dh = c/H is the Hubble distance. 
Specializing to a flat Universe in Fig. 2 we display m(z) in Eq. (5.18) as a function of the 
redshift z for the three models. The recent SN la sample provides measurements of the 
luminosity distance out to redshift z ~ 1.7 (SN 1997ff). In Figure 2 we also shows the 
Hubble diagram of effective rest-frame B magnitude corrected for the width-luminosity 
relation, as a function of redshift for the 42 Supernova Cosmology Project high-redshift 
supernovae [22], along with the 18 Calan/Tololo low- redshift supernovae [2B]. We also 
display the 10 High-Z Supernova search Team supernovae with 0.16 < z < 0.62 [21], and 
the recent SN 1997ff at z ~ 1.7 [2H]. In the case of the farthest known supernova SN 
1997ff we also report the revised value of effective rest-frame B magnitude corrected after 
correction for lensing |26j . 

Even though we do not attempt to best fit the supernova data, from Fig. 2 it is evident 
that the model with Q mQ = 0.3, f^Ao = 0.0 , Q w = 0.7 is excluded. On the other hand, 
both the standard model and the model with Q mQ = 0.3, f^Ao = 0.6 , fl w = 0.1 are 
consistent with the high-z supernova data. " This can be better appreciated from Fig. 3 
where we display the difference between data and the fiducial standard model fl m0 = 0.3, 
^ao — 0.7, Q w = 0.0. Thus we see that recent high-z supernova data allow that a 
small, but cosmological important, part of our Universe consists of an almost ideal gas of 
ferromagnetic domain walls. 



6. CMBR constraints 



As discussed in Section I, in general we expect that a network of domain wall extend- 
ing over cosmological distance could lead to severe distortions of the Cosmic Microwave 
Background Radiation. Thus, it is important to investigate the influence of the cosmo- 
logical gas of massive ferromagnetic domain walls on the Cosmic Microwave Background 
Radiation. If domain walls are present we have an additional temperature anisotropy: 

= A$ , (6.1) 



Hit is worthwhile to stress that the perfect gas approximation (which implies that a network of FDWs 
behaves as a radiation gas) cannot be valid when interactions with the plasma are taken into account. 
Indeed, if the equation of state were valid at all times then, due to the large value of relativistic matter 
(Cl w = 0.1), the equilibrium between matter and radiation would take place later respect to the cosmo- 
logical standard model and the structure formation would proceed at too small z. It is clear that this 
would be in contrast with the standard analysis of structure formation (see e.g. On the other hand, 

assuming an equation of state more stiff for very high redshifts, such high values for tt w are not a priori 
excluded. 
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where $(r) is the Newtonian potential of the wall network. To evaluate the temperature 
anisotropy Eq. (6.1) we follow the analysis performed in Ref. [27]. We find 

^ = *{(L))n», (6.2) 

where (L) is the average distance between walls and n is the number of walls per horizon. 
As concern the exponent u, it turns out that this exponent depends on the network walls 
configuration. According to Ref. [27J we consider values of v between v — 1 (regular 
network) and v — 3/2 (non-evolving network). In our case, the Newtonian gravitational 
potential for a FDW in the (x,y) plane is <fi(z) = 4ttGn\z\<j w **. Moreover we have 
n = d H / (L), so that from Eq. ()4.2|) we obtain: 

^r ) =^G N o*-^-") d H Pr 1 • (6-3) 

The quantities in Eq. (6.3) need to be evaluated at the present time to- Let us observe 
that: 

, / x / s f t0 dt f 1 dr . . 

d H (to) = R(to) / = h / , 6.4 

Jo R{t) Jo t{t) 

where r(r) = R(t)/R{t ), r = t/t , and tf 



>3 



cr w {ta) = a w (t ew ) f ^/^V ) = M*<™) r'L , (6.5) 



t ew denoting the time at the electroweak phase transition. 
Now, we recall that r(r) is solution of the equation: 



d ^ = t H [tlmor- 1 + n w0 r~ 2 + Q A0 r 2 ] 1/2 . (6.6) 



Putting f(r) = (n m0 r 1 + Q w qt 2 + Vt^r 2 ) 1 / 2 we have 



So that we can write: 

ST 



(T) 

where 



a u [G N a w (t ew )t 3 J w 2 H 1 Q /2 } 2 -', (6.8) 



**Given the stress tensor Tjf = diag(p,—pi,—p2,—p3), the Newtonian limit of Poisson's equation is 
V 2 = 4ttGn(p + Pi + P2 + P3), where <fi is the Newtonian gravitational potential. For a thin FDW 
p = 8(z)a w , p\ = p-2 = p/2 and pa = 0. Thus, V 2 </> = 8itGnS(z)<7 w . 

^From energy conservation we have that p w ~ i?~ 4 ; since p w ~ a w /(L) and (L) ~ i?, we get <r w ~ i? -3 . 
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Eqs. (|6.8I6.9|) give an approximate estimate of ST/(T). 

Taking G N ~ (10 19 GeV)" 2 , a w {t ew ) ~ 10 6 GeV 3 , t ew ~ 10 13 GeV -1 , H ~ 10~ 42 GeV|T7], 
and assuming Q m = 0.3, fl w = 0.1, = 0.6, we obtain: 

^ ~ 10- 33 for v = l, (6.10) 

— ~ lO- 1 ^" 1 for u = 3/2. 

Even in the worst case v = 3/2 the contributions to the large-scale anisotropy of the 
microwave background radiation is completely negligible given the observed value ST/T ~ 

io" 6 Hang. 



7. Conclusions 



Let us conclude by stressing the main results of the paper. From the constraint that 
domain walls must not destroy the observed isotropy of the Universe, it is widely believed 
that domain walls that might have been produced in the standard particle physics phase 
transitions in the very early Universe have been eliminated by some mechanism. We have 
discussed a new mechanism which allows domain walls produced during the primordial 
electroweak phase transition. Indeed, the effective surface tension of these domain walls 
can be made vanishingly small due to a peculiar magnetic condensation induced by fermion 
zero modes localized on the wall. As a consequence, the domain wall acquires a non 
zero magnetic field perpendicular to the wall, and it becomes almost invisible as far as 
gravitational effects are concerned. We find that in the perfect gas approximation the 
domain wall network behaves like a radiation gas. The analysis of the recent high-redshift 
supernova data suggests that a small component of our Universe could be composed of 
this gas of ferromagnetic domain walls. Finally, we showed that ferromagnetic domain 
wall gas induces a completely negligible contribution to the large-scale anisotropy of the 
microwave background radiation. 

Undoubtedly, the presence of a network of massive FDWs can have influenced physical 
processes that took place in the early Universe. It should keep in mind that in studying 
effects occurred in the early Universe (say before or during the formation of large scale 
structures) the perfect gas approximation in no longer valid and a more complicated 
analysis of the interaction of walls with the primordial plasma must be in order. We 
reserve such subject matter to future investigations. 
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A Appendix 



We are interested in the Dirac equation with a kink wall and in presence of the electro- 
magnetic field Ap(x): 

(ift-e4-gy<f>)il) = 0. (A.l) 
To solve Eq. (A.l) we assume that 



if, = (ifi - e 4 + 9y 4>) 
Inserting Eq.(A.2) into Eq. (A.l) we readily obtain: 

_ 2ie A p + e 2 A 2 _ ^2 g 2 j rii p g _ t^F ) $ = 



(A.2) 



(A.3) 



with £ = gyv, g(z) = tanh(z/A). It is easy to see that the solution of Eq. (A.3) factorize 
as: 

$(x,y,z,t) = f{x,y)x(z,t), (A.4) 

where / is a scalar function which describes the motion transverse to the wall and \ 
is the spinorial part of the wave function of fermions localized on the wall. Putting 

x(z,t) = e~ lEt x{z) we get for / and x : 



(d 2 + d 2 - 2ieBxd y - e 2 B 2 x 2 + (3) f(x, y) = 



(A.5) 



(d 2 z + i^ 3 d z g - ^ 2 9 2 + 2eB 7 V + a) X (z) = 0, 
where a and (3 are constants subject to following relation: 

E 2 = a + (3. 

The solution of Eq. (A.5) is: 

f n ^y) = A n e-^ 2 + ^H n (C), 

where 

y/2 n n\ V eBJ 

H n (x) are Hermite polynomials and 

P= |e|B(2n + l), n = 0,l,2,... (A. 10) 

In order to solve Eq. (A. 6), we expand x( z ) i n terms of spinors w± 2 eingenstates of 7 3 

X(z) = <p 1 ± 2 u 1 /. (All) 

Using the standard representation for the Dirac matrices [2H] we find: 



(A.6) 
(A.7) 

(A.8) 
(A.9) 



u: 



( 1 \ 



±i 

V o J 



u_ 



( \ 
1 


V Ti ) 



(A.12) 
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It is straightforward to check that: 



1,2 1,2 / a 1 o\ 

7 «± = m t ' ( A - 13 ) 

7 1 u| = =R^, 

2 1,2 _,_ 2,1 

T«± = =P^ , 

3 1,2 ,. 1,2 

7 = ira^ . 



Taking into account Eqs. (A. 6), (A. 11) and (A. 13) we obtain: 

{d 2 z =F id z g - i 2 g 2 + eB + a) <f>l(z) = 0, (A.14) 



[d\ =F ^ - £ 2 <? 2 -eB + a) <f> 2 ± (z) = 0. (A.15) 
It easy to see that there are localized states if 

a = -eB for (A.16) 
a = +eB for 0|. (A.17) 

Inserting Eqs. (A. 10), (A.16), (A.17) into Eq. (A. 7) we get the energy spectrum for local- 
ized states: 

E n = ^2n\e\B n = l,2,... (A.18) 
We can now solve Eqs. (A.14) and (A.15) whit the constrains (A.16) and (A.17); we find: 

cj>¥(z) = N exp (±£ J*dz>g(z')^, (A.19) 

where the normalization constant iV is evaluated below. The solutions 0+ 2 do not corre- 
spond to localized states and will be discarded. Inserting Eqs. (A.19), (A. 11), (A. 8), and 
(A. 4) into Eq. (A. 2) it easy to recover our Eqs. (3.2)-(3.5). Finally, the normalization 
constant N can be obtained by imposing the normalization conditions: 



d 3x Vv!, P , s ( x ) Vvy,s'( x ) = Kn< S ss 'S(p - p') . (A.20) 
It is straightforward to obtain: 

2 = y^\B 

8ttAB(£A,1/2)£ 2 ' 1 ' ; 
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Figure 1: Evolution of expansion parameter for three different models. 
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Figure 2: Hubble diagram for Type la supernovae data. The meaning of symbols is as 
follows. Calan/Tololo supernovae (full triangles down), Supernova Cosmology Project 
high-redshift supernovae (full circles), High-Z Supernova search Team supernovae ( full 
triangles up), SN 1997ff (full square) and SN 1997ff after correction for lensing (open 
square). Full line corresponds to the standard model fl m = 0.3, £7a = 0.7 and Q w = 0, 
dashed line to Q m o = 0.3, ^ao = 0.0 , Q w = 0.7, dotted line to f2 m0 = 0.3, Q\o = 0.6 , 
tt w = 0.1. 
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Figure 3: The magnitude residuals between data and models f2 m o = 0.3, ^ao = 0.0 , 
Q w = 0.7 (dashed line), f2 m0 = 0.3, f^Ao = 0.6 , Q w = 0.1 (dotted line) from the standard 
model Q m = 0.3, ^a = 0.7 and Q w = (full line). 
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